Draft version September 28, 2009 

Preprint typeset using I^'T^]X style cmulatcapj v. 08/22/09 



DUST PRODUCTION AND MASS LOSS IN THE GALACTIC GLOBULAR CLUSTER NGC 362 

Martha L. Boyer\ Iain McDonald^, Jacco Th. van Loon^, Karl D. Gordon^ Brian Babler^, Miwa Block'', 
Steve Bracker"', Charles Engelbracht"*, Joe Hora\ Remy Indebetouw'', Marilyn Meade'', Margaret Meixner\ 
Karl Misselt^, Joana M. Oliveira^, Marta Sewilo\ Bernie Shiao\ and Barbara Whitney' 

(Received; Revised; Accepted) 
Draft version September 28, 2009 



o 
o 

(N 

D 

00 
(N 

00 

6 



> 
in 

0^ 
O 
On 
O 



ABSTRACT 

We investigate dust production and stellar mass loss in the Galactic globular cluster NGC 362. Due 
to its close proximity to the Small Magellanic Cloud (SMC), NGC 362 was imaged with the IRAC 
and MIPS cameras onboard the Spitzer Space Telescope as part of the Surveying the Agents of Galaxy 
Evolution (SAGE-SMC) Spitzer Legacy program. We detect several cluster members near the tip of 
the Red Giant Branch that exhibit infrared excesses indicative of circumstellar dust and find that 
dust is not present in measurable quantities in stars below the tip of the Red Giant Branch. We 
modeled the spectral energy distribution (SED) of the stars with the strongest IR excess and find a 
total cluster dust mass-loss rate of 3.0^^ 2 x 10~^Mq yr"-'^, corresponding to a gas mass-loss rate of 

8.6l3'4 X XQ^^Mq yr^'^, assuming [Fe/H] = —1.16. This mass loss is in addition to any dust-less mass 
loss that is certainly occurring within the cluster. The two most extreme stars, variables V2 and V16, 
contribute up to 45% of the total cluster dust-traced mass loss. The SEDs of the more moderate stars 
indicate the presence of silicate dust, as expected for low-mass, low-metallicity stars. Surprisingly, 
the SED shapes of the stars with the strongest mass-loss rates appear to require the presence of 
amorphous carbon dust, possibly in combination with silicate dust, despite their oxygen-rich nature. 
These results corroborate our previous findings in lo Centauri. 

Subject headings: globular clusters: individual (NGC 362) — stars; mass loss — circumstellar matter 
— stars:winds, outflows — infrared: stars — stars: AGB and post-AGB 



1. INTRODUCTION 

Stellar mass loss and dust production remain two of 
the most critical, yet least understood, aspects of stel- 
lar evolution models, especially in low-mass population 
II stars. High-mass Asymptotic Giant Branch (AGB) 
stars (up to 8 Mq) evolve and lose mass on very short 
timescales, resulting in short bursts of dust input into 
the Interstellar Medium (ISM) following episodes of star 
formation. Low-mass stars are more numerous than their 
high-mass counterparts and live much longer, resulting in 
a more sustained dust input into the ISM. Constraining 
these aspects is crucial for determining the impact of low- 
mass stars - the most numerous stars in the Universe - 
on galaxy evolution. Observat ions wi th the Spitzer Space 
Telescope ( Werner et al. 200 l:lGehrz et al. 2007) and the 
AKARI telesco'pe (|Murakami et al.ll2007,l have produced 
many studies of infrared photometry of resolved stellar 
populations in globular clusters (GCs). Despite these 
new studies, the sample of observations is still quite 
small, consisting of only a handful of GCs and creating 
almost as many new questions as it has answered. 
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lOrigha et al.l ()2007n argued on the basis of their analy- 
sis of Spitzer data for 47 Tuc that dust-accompanied mass 
loss occurs along the entire Red Giant Branch (RGB), 
down to the Horizontal Branch (HB). However, this 
was not confirmed with AKARI observations (|Ita et all 
[2OOI . In other GCs, dusty mass loss is only seen near the 
very tip of the RGB (jBover et al.ll200a iMcDonald et all 

Eoog). 

Intracluster dust clouds should form as a result of 
dusty mass loss from many stars during the time between 
Galactic plane crossings. Searches for this dust in a large 
sample of GCs using Spitzer and AKAR I have resulted 
in n o successful/significant detections (jBarmbv et al.l 
120091: iM atsunaga et al. 2008), with the exception of M15 
(jBover et al.l I200Q) . This lack of intracluster medium 
(ICM) dust in most clusters provides a mystery as to the 
fate of dust produced by evolved stars in GCs. 

1.1. NGC 362 

We present an infrared study of NGC 362, a bright GC 
in the southern sky. As one half of the classical "second- 
parameter" GC pair with NGC 2 88, NGC 362 is a rel- 
atively well-st udied cluster (e.g.. iBellazzini et al.l 120011 
iCatelan et al.l |2001). The wide availability of multi- 
wavelength and astrometric data, along with the cluster's 
large coeval stellar population, a metallicity intermedi- 
ate to 47 Tuc and the bulk populati on of w Cen ([Fe/H] 
= -1.16 for NGC 362. lHarrisl ll996h^ and proximity to 
us (8.5 kpc. IHarrisl[T996i r make it an ideal candidate for 
a study of dusty mass loss on the AGB and upper RGB. 

^ The lHarr is (1996) catalog was update in Februrary 2003. See 

|http: / /www.physwww . physics . mcmaster .ca/'^harris/ mwgc . dat ^ 



2 



Boyer et al. 




NGC 362; 3.6 micron (blue), Q micron (green), 24 micron (red) 



TABLE 1 
NGC 362 Properties 



Parameter 



Value 



Source 



Right Ascension (J2000) 0l'^03™14!27 

Declination (J2000) -70°50'53^'6 

Distance (kpc) 8.5 3 

[Fe/H] -1.16 3 

M (Mq) 3.78 X 10^ 2 

Rcorc (') 0.17 3 

-Rhalf mass (') 0.81 3 

E(B - V) (mag) 0.05 3 

(m - M)v (mag) 14.65 3 

Tc (yr)'' 3 X lO'^ 4 

iVHB 1.6 X 102 1 

VeBcO (kms-i)!^ 46.7 5 

My (mag) -8.35 3 

Lboi (Lq) 4.54 X 10^ 3 

Heliocentric radial velocity (kms~^) 223.5 it 0.5 3 



Note. 



(11996 



Sources: (1 ) 

iGnedi n fc Os trikeH 
(4) [Odcnkirc hen et al 



iBarmbv et al.l 
1997D. (3) 
11993), (5) 



McLaughlin fc van der Mareli (^2005!) 

Tc is the time since the last Galactic plane-crossing. 
^ Wcsc,o is the escape velocity at i? = 0. 

Gompared to NGC 288, NGG 362 has a red HE, de- 
spite having a very similar metallicity to the former 
(i.e., the "first parameter"). There are several theories 
for this discrepancy, including (and possibly combining) 
cluster age, helium abundances, and RGB mass loss, 
which may be affected by the central co ncentration of 
the c luster or other environmental factors (| Cat elan et al.l 
i2001l ). A recent study also suggests environmental condi- 
tions during formation as a second-parameter candidate 
(jFraix-Burnet ct al. 2009). 

NGC 362 appears to be devoid of an IGM. 
iBarmbv et aH p009. ) find an upper limit of 6.3 x IQ-^ M© 
of dust, more than two orders of magnitude less than 
predicted based on the number of HB stars and the 
time since the last Ga lactic plane-crossi n g (tv , see 
Table ffl). In a dditio n. iCrindlav fc Lilled |l977i ) and 
iHesser fc Shawli (Il977| ) searched for and found no ion- 
ized intracluster medium. Despite the apparent lack of 
material collected from mass-losing stars in the clus- 
ter, two studies have identified a small population of 
sources with infrared excesses attributed to circumstellar 
dust and mass l oss (Origlia ct al. 2002; Ita ct al. 2007), 
and iMcDonald fc van LoonI (j2007l) have identified clus- 
ter stars with Ha profiles that indicate mass loss in the 
absence of dust. 

2. OBSERVATIONS AND DATA REDUCTION 

A 3-color image composed of 3.6, 8, and 24 /im In- 
frared Array Camera (IRAC) and Multiband Imaging 
Photometer (MIPS) images of NGC 362 is presented in 
Figure[T] The cluster's near-juxtaposition with the Small 
Magellanic Cloud (SMC) resulted in serendipitous obser- 
vations of NGC 362 with Spitzer as part of the Survey- 
ing the Agents of Galaxy Evolution Spitzer Legacy Pro- 
gram (SAGE-SMC; K.D.Gordon 2009, in preparation). 
NGC 362 lies approximately 2° north from the center of 
the SMC bar, placing it near the edge of our 5° x 5° 
Spiiteer observations. The cluster was covered at 3.6, 4.5, 
5.8, 8, 24, and 70 fj,m to well outside of its half-mass ra- 



FlG. 1.— 3-color IRAC and MIPS image of NGC 362. Blue is 
3.6 fim, green is 8 fim, and red is 24 fiva. The northern most 
portion of the SMC bar lies approximately 1° to the south. 



dius of 0.81' (IHarrislll996n . SAGE-SMC 160 fim images 
did not include the cluster. All data and analysis pre- 
sented here are confined to within 6' of the cluster center 
to minimize contamination from SMC sources. 

Observations consist of two epochs of IRAC images 
separated by 3 months and one epoch of MIPS images. 
The 1 a sensitivities are 2.3, 3.4, 19.1, 20.5, and 28.1 /xJy 
for 3.6, 4.5, 5.8, 8, and 24 ^m, respectively. Short expo- 
sure times (12 s) ensure that saturation of cluster stars is 
not an issue. Angular resolutions for IRAC wavelengths 
range from 1.7" at 3.6 /im to 1.9" at 8 /im and increase 
to 5.8" for MIPS 24 /im. For more details regarding data 
acquisition and reduction for the SAGE-SMC program, 
see K.D.Gordon, et al. (2009, in preparation). 

Photometry should be reasonably complete to well be- 
yond the HB (MjjB ~ — 1 mag at 3.6 /im), as demon- 
strated by a steady increase in source counts to Afa.gpm ~ 
2 mag in Figure [21 While photometric completeness 
tests have not yet been performed on the SAGE-SMC 
data, the completeness limits beyond the crowded in- 
ner 1' of the cluster are likely similar to those for the 
SAGE Large Magellanic Cloud (SAGE-LMC) IRAC data 
(pvlcixncr et_al. .201)6) since the observations were simi- 
larly designed. At 3.6 /im, the average photometric com- 
pleteness limit for SAGE-LMC is f»16 mag (Ms.g^m ~ 
1.4 mag for NGC 362). The luminosity function in the 
inner 1' of the cluster is truncated just below the HB 
(Fig.[2|); this is a result of incompleteness due to crowd- 
ing. The fact that the luminosity function at brighter 
levels does not differ in the cluster core from the field 
means that crowding is not important for stars brighter 
than the HB. 

Only three point sources are detected at 70 //m, and 
only one of these (designated here as si 1) is also detected 
in IRAC. This source is discussed further in Section [3T2l 

2.1. Ancillary Data 

To aid in estimating the temperatures and luminosi- 
ties of the cluster stars, we collected literature photome- 
try spanning the optical and the near-IR. NGC 362 was 
observed in the Magellanic C louds Photometric Survey 
(MCPS; Zarit skv et all [200^ the Deep Near-Inf rared 
Southern Sky Survey (DENIS: lE^chtein et al.lil997l ). the 
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Naval Observatory Merged Astrometric Dataset (NO- 
MAD; iZacharias et all l2005f) and the Tw o Micron All 
Sky Survev (2MASS: I Skrutskie et al.l[200l) . 

The MCPS catalog includes UBVI photometry, but 
does not cover the inner wO.B' of the cluster core, where 
four candidate mass-losing stars identified with Spitzer 
data reside (see Section ^ . We obtained optical pho- 
tometry of stars in the cluster core from the DENIS and 
Nomad catalogs, which include /-, J-, and ifg-band pho- 
tometry and V-, and i?-band photometry, respec- 
tively. 2MASS includes JHKg photometry and probes to 
the cluster center. This collection of optical and near- 
IR data complements the 3.6 /zm Spitzer data in both 
completeness and angular resolution, with 93% of the 
3.6 /im sources matched in at least one of these four an- 
cillary catalogs. The remaining 7% of 3.6 /zm sources 
without optical or near-IR counterparts have a mean 
3.6 /im magnitude of 16.5 mag, with a standard devi- 
ation of 1.3 mag, or three magnitudes fainter than the 
HB. All of the sources investigated in detail in this study 
(Section[3]) are detected in at least one near-IR or optical 
catalog. 

In addition. Ilta et ahl ()2007l ) obtained infrared photom- 
etry at 2.4, 3.2, 4.1, 7.0, 9.0, 11.0 15.0, 18.0, and 24 ^m 
with the Infrared Camera (IRC: lOnaka et aI1l2007f) on- 
board the AKARI telescope (see Fig. 3 from Ilta et al.l 
120071 photometry obtained through private communica- 
tion). IRC images have angular resolutions ranging from 
0.9" at 2.4 ^Tci to 9" at 24 /im, making stellar blending 
a problem in the most crowded regions of the cluster. 

NGC 362 was al so observed at 24 and 70 /im by 
IBarmbv et aP ((20091) . While none of the candidate mass- 
losing stars identified in Section [3] are detected at 70 /im, 
the supplementary 24 /im data are useful for estimating 
dust compositions and mass-loss rates of individual stars 
(see Section [33]) . 

A Hubble ACS image was obtained with permission 
from the ACS Survey o f Galactic globular clusters team 
(|Saraiedini et al.ll2007f) . The high resolution of this im- 
age was helpful in determining when stellar blending in 
Spitzer and AKARI images is potentially a problem. 

2.2. Cluster Membership 

Proper moti ons in NG C 362 were o btained by 

IZacharias etldl (|2005[ ) and iTucholkd (fl99l . The lat- 
ter computed a membership likelihood by fitting the dis- 
tribution of stars in /i^cosJ, /i^ space with a sum of 
two bivariate Gaussians representing NGC 362 and field 
stars. We use these probability estima tes to eliminate 
non-m embers. For proper motions from IZacharias et al.l 
(|2005f ), we employ a conservative cut-off of >35 mas yr~^ 
to eliminate probable cluster non-members. 

The heliocentric radial velocity (i^rad) of NGC 362 is 
223.5 ± 0.5 km s'^ (|Harrislll996n . iFischer et all ([l99l ) 
measured radial velocity for 215 stars around the center 
of the cluster, and confirmed t wo stars in their sample as 
radial velocity non-members. iHarris fc Zaritskvl ()2006l ) 
measure the radial velocities of red giants in the SMC as 
146 ± 28 km s"\ and find that it is rare for an SMC 
star to have a radial velocity larger than 200 km s~^. 
However, they did not measure stars at positions near 
NGC 362, and we therefore note that it is possible, al- 
though rare, for SMC stars to have radial velocities sim- 




-8 -6 -4 -2 2 4 

Fig. 2.— Luminosity function of NGC 362 (within R = &), in- 
cluding both epochs of IRAC data. The sohd and dotted Hnes illus- 
trate 0.2 and 0.5 magnitude bins, respectively. The tip of the Red 
Giant Branch (TRGB) is located near Ma.e^m ~ — 6.2 mag, and 
the Horizontal Branch (HB) is located near Afa.e^m ~ — 1.0 mag. 
Photometry should be reasonably complete to approximately 3 
magnitudes fainter than the HB at 3.6 /im. In the inner 1' of 
the cluster (gray line) crowding limits source completeness to mag- 
nitudes brighter than the HB. 

ilar to NGC 362 due to the rotation of the SMC (for 
relatively high-mass stars) and its velocity dispersion (in 
particular for low-mass stars). 

Membership information is available for 68% of the 
sources detected by Spitzer. We suspect that many of 
the sources without membership information (especially 
those with extreme IR excess) are in fact background 
galaxies that are not detected in the optical (see Sec- 
tioning). 

3. MASS-LOSING STARS IN NGC 362 
3.1. Previously Identified Mass-Losing Stars 

Three so urces in NGC 362 (designat ed as xOl, x02, 
and x03 by iMcDonald fc van Loonll2007f) were ide ntified 
as having infrared excess by lOriglia etall ((200l usmg 
12-/im Infrared Space Observatory (ISO) data. The co- 
ordinates of source xOl coincide with a Spitzer source 
that shows 24 /xm ex c ess ( source s04, see Table [5]). 
IMcDonald fc van LoonI (|2007l ) showed xOl to have a red 
Ha line emission wing and blue-shifted line absorption 
core, indicative of an outflow with a mass-loss rate of 
'-~^10~^ Mq yr~^. Two other IR excessive point sources 
(sources s03 and s05) lie immediately to either side of this 
source, and while resolved in Spitzer, all three sources are 
likely blended together in ISO images. 

ISO source x02 corresponds to three different objects, 
two of which (x02a and x02b) are detected with Spitzer 
and show little to no IR excess despit e having blue- 
shifte d Ha absorption cores (McDona ld fc van LoonI 
[2007). The third source (x02c) is unresolved by Spitzer. 
Although apparently exhibiting a lack of IR excess, we 
note that x02a and x02b are located in the most dense 
region of the cluster where source confusion in Spitzer is 
severe and reliable photometry is difficult. 

ISO source x03 is not resolved at 24 /im with MIPS, 
but its IRAC colors do indicate a possible (slight) 8 /im 
excess. The AKARI observation of x03 measures 3.2, 
4.1, 7, and 9 /im fluxes that are four times brighter than 
the Spitzer 3.6, 4.5, and 8 fim fluxes, likely due to stel- 
lar blending in AKARI since x03 is located in a dense 
region of the cluster. Despite this discrepancy, AKARI 
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fluxes still show x03 to have a very slight IR excess. A 
strong red emission wing in the star's Ha line and blue- 
shifted absorption c ore yield a potential mass-loss rate 
of -lO-^ Mq yr-i (|McDonald fc van Lo on 2007). The 
lack of reliable photometry redward of 8 /im prevents us 
from measuring a mass- loss rate in Section 13.51 

In addition to the three ISO sources, 
iMcDonald fc van LoonI pOOl obtained VLT/UVES 
spectra of 10 other stars in NGC 362, all of which have 
spectra typical for oxygen-rich red giant stars. Five of 
these sources show IR excess at 24 /im (see Table [2] 
and Section 13. 5|) . Fitting of the Ha profiles suggests 
mass-loss rates ranging from 10~^ — 10~^ Mq yr~^. 
Source oOl [Spitzer source s0 6) has strong I R excess, is 
variable (V16, P = 138 davs: lClementl ll997H . and shows 
strong molecular bands and emission in Ha. 

The remaining fiy e sources from 

IMcDonald fc van LoonI ()2007D are also detected in 
IRAC. None of these sources show significant 8 or 24 /^m 
excess, the exception being o09, which shows slight 
8 /^m excess ([3.6] — [8] — 0.18, Ma.g^m = —5.1 mag). 
It is interesting to no t e that all sources observed by 
IMcDonald fc van LoonI (|2007D that show IR excess m 
Spitzer or ISO data also show red Ha line emission 
wings. The four stars (o02, o03, o04, and olO) that 
have no IR excess show only b l ue Ha line emission 
wings. IMcDonald fc van LoonI (j2007f ) suggest that 
stars showing only blue emission wings have heated 
chromospheric material restricted to radii less than 
2 Rq , resulting in the red emission being blocked by the 
stellar disk. In this scenario, the chromospheres of stars 
showing significant IR excess will be extended to large 
e nough radii t o see a red emission in Ha. 

iSmith et al.l (119991 ) obtained an optical spectrum of 
another NGC 362 variable from iClement (1997), V2 
[P = 90 d ays). V2 is Spitzer somce s02 (Tabled Sec- 
tion HSH]), which has strong IR excess. Strong Ha and 
H/3 emission in V2 indicate chromospheric activity, and 
a strong Li I A6707 feature indicates a lithium overabun- 
dance. Low-mass giant stars have typically destroyed 
and diluted their photospheric lithium (which is brought 
to the surface via the first dredge- up), so those showing 
a lithium overabundance are rare, although not unheard 
of. 

Eight sources in NGC 362 showing strong IR ex- 
cess (F24^m/F7^m > 1, whcre zero-e xcess star s have 
^2A^^ln/^^^^m ~ 0.1 ) wcrc identified bv llta et all (|2007D 
using AKARI data. These eight sources are located in 
the region of the IR color-magnitude diagram (CMD) 
that is now attributed to b ackgr ound galaxies (e.g., 
Blum et al. 2006; Bover et al.l l2008L 120091 : iBolatto et all 
2007: McDonald ct al. 200^ although there may also 
be a small amount of contamination from background 
evolved stars in this region. Several of the cluster mem- 
be rs showing IR ex cess in Spitzer data are also detected 
bv llta et al.l ()2007l ). but these sources are not the focus 
of that particular study. 

3.2. Identifying Mass-Losing Stars with Spitzer 

A CMD showing Spitzer [24] versus [8] - [24] is pre- 
sented in Figure HI Confirmed non-members are repre- 
sented with black dots, and sources without member- 
ship information are marked with gray dots, the major- 
ity of which are likely background galaxies. Stars on the 
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[8] - [24] 

Fig. 3. — [24] versus [8] — [24] color-magnitude diagram. Lines 
connect points from each IRAC epoch. Five sources are identified 
as candidates of strong dusty mass-loss, and five others are candi- 
dates for moderate dusty mass loss. One of the moderately dusty 
candidates is a likely cluster non-member (source sOI). 

RGB have average colors of [8] — [24] — mag, while 
stars near the bright end of the RGB begin to shift to 
[8] - [24] > mag. 

Table [2] lists ten sources with [8] — [24] > mag lo- 
cated in the region of the CMD above or near the tip of 
the Red Giant Branch (TRGB; Ms.e^m ~ - 6.2 mag, 
based on the drop-off in the lumino sity funct i on, s hown 
in Figure H]) . Proper motion from iTucholkd (|1992i ) in- 
dicates that source sOl is a likely non-member, but the 
remaining nine sources are confirmed radial velocity or 
proper motion cluster members. Five of the ten sources 
(s02, s05, s06, s07, and s08) are either very bright or 
have strong 24 /zm excess, and we designate them as can- 
didate strong mass-losing stars. The four more moder- 
ate confirmed member stars are potentially forming dust 
in smaller quantities and possibly also losing mass at a 
moderate rate. 

In addition to the ten sources listed in Table |2l source 
sll (R.A. = 1'>02"M8!66, Dec. = -70°45'22'.'6) is very 
red and bright, with [8] - [24] > 3.5 and M24^m < - 
7 mag. This source is the only IRAC source also detected 
at 70 /im, and it is most likely an unresolved background 
galaxy, despite being more than a magnitude brighter 
than the general locus of other unresolved background 
galaxies. Its position in the [3.6] versus [3.6] — [8] CMD 
(Fig. H]) and its inferred bolometric luminosity if it were 
a cluster member (see Section 13.3! and Figure [7]) place it 
among the other background galaxies, far from the mass- 
losing AGB stars. Moreover, this source's SED (Fig. [5]) 
rises strongly beyond 3 /im, re sembling those o f galax- 
ies in the SINGS sample (e.g.. lDale et al.ll2006l ) and of 
other background galaxies identified in Figure |31 While 
we cannot confirm this source's status as a background 
galaxy without a spectrum or high-resolution image (Sll 
falls outside of the ACS field of view), we proceed under 
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TABLE 2 

Candidate Strong and Moderate Mass-losing Stars in NGC 362 



Source R.A. (J2000) Dec. (J2000) Ms^m M24pm T L Notes and 

ID (h m s) (° ' ") (mag) (mag) (K) (Lq) Alternate Designations^ 



Moderate Mass-Loss Candidates 



sOl 


01 03 35.75 


-70 50 52.4 


-6.17±0.03 


-6.38±0.02 


3950 


1656 


non-member 


s03 


01 03 20.07 


-70 50 55.0 


-6.09±0.03 


-6.30±0.02 


4339 


2184 


o07, post-AGB? 


s04 


01 03 19.07 


-70 50 51.1 


-6.07±0.03 


-6.07±0.03 


3823 


1551 


xOl 


s09 


01 03 10.67 


-70 50 54.4 


-6.08±0.03 


-6.18±0.03 


4226 


1932 


o08, post-AGB? 


slO 


01 02 43.34 


-70 48 47.8 


-5.81±0.03 


-5.98±0.03 


3975 


1363 





Strong Mass-Loss Candidates 



s02 


01 03 21.85 


-70 54 20.2 


-6.73±0.03 


-7.15±0.02 


3907 


1826 


V2 (P = 90 days) 


s05 


01 03 17.29 


-70 50 49.1 


-6.26±0.05 


-6.68±0.02 


4058 


2280 


o05a/o05b 


s06 


01 03 15.08 


-70 50 31.5 


-7.28±0.03 


-7.45±0.02 


3962 


3106 


oOl, V16 (P = 138 days) 


s07 


01 03 13.62 


-70 50 36.5 


-6.40±0.03 


-7.20±0.02 


3343 


1402 


o06 


s08 


01 03 12.60 


-70 51 00.9 


-5.31±0.02 


-6.08±0.02 


3682 


1842 


stellar blond? 



Note. — The prefix of the source IDs is short for "Spitzer'\ Sources are numbered in order of decreasing R.A. 
Designations fro mlMc Donald fc van LoonI II2007I ) have a prefix "o" or "x", followed by a number. Source s04 (xOl) 
was also detected by [OrTelT^'et'^.l 120021 ). Photosph ere temperatures derived here for s03 and s09 are warm, possibly 
indicating that these are post-AGB stars. However, IMcDonald fc van LoonI 120071) estimate temperatures that place 
both stars firmly on the RGB. 
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Fig. 5. — SED of source sll, a suspected background galaxy. 
Black pluses, blue circles, and red circles represent optical/Near- 
IR, Spitzer, and AKARI data, respectively. Wavelength error bars 
correspond to the filter widths. The scatter in the optical and 
near-IR, photometry is likely due to mis-matching of sources, as 
foreground optical sources may lie directly along the sightline to 
the galaxy. Four other sources from the galaxy locus in Figure [3] 
are also shown (green diamonds, purple triangles, orange squares, 
and blue asterisks), plotted with offset fiuxes for visual clarity. The 
SED of sll resembles that of the other four background galaxies. 



Fig. 4. — [3.6] versus [3.6] — [8] color-magnitude diagram. Lines 
connect points from each IRAC epoch. The sources identified as 
candidate stars with strong and moderate mass loss in Figure |3] arc 
marked with red circles and purple diamonds, respectively. The 
crowded region towards the TRGB is inset at the top right. The 
source at [8] — [24] > 3.5 and A/24pm < — 7 mag in Figure |3] (source 
sll) is located far from mass- losing AGB stars and near unresolved 
background galaxies in this CMD. Source sl2 may be a background 
SMC carbon star. 

the assumption that it is not a member of NGC 362 and 
exclude it from our analysis. 

It is clear from the Spitzer data that dusty mass loss in 
NGC 362 is confined to the upper RGB/ AGB. There are 
a handful of potential member stars with moderate 8 fxin 
excess at —2 > Ma.g/jm ^ —4 mag visible in Figured 
All of these sources are very near the cluster center, and 
upon visual inspection, all but one are clearly blended 
with other sources in all of the Spitzer data, potentially 
creating artificially red [3.6] — [8] colors if the flux is not 



accurately extracted at one or both wavelengths. Indeed 
it is likely that the flux of these sources is inaccurate 
since these are only moderately bright stars located in 
an extremely crowded region of the cluster. Given that 
the mean uncertainty in their [3.6] — [8] colors is 0.17 mag 
(0.2 > [3.6] — [8] < 0.3 mag), these sources could easily 
be members of the non-dusty RGB. 

The only potentially unblended source in this re- 
gion of the CIVID (sl2; R.A. = l'^03'"ll'!59, Dec. = 
-70°50'23'.'5) is located at [3.6]_- [8] « 0.47 ±0.14 and 
A^3.6Aim ~ — 4 mag in Figure [H T his source was desig- 
nated a radial velocity member by iFischer et al.l ()1993D 
('^rad = 224.8 km s^ ), but it has no known proper mo- 
tion measurements. IR data for this source are somewhat 
ambiguous (Fig. [5]), and it is very likely that either the 
IRAC photometry is inaccurate or this source is actually 
a background star belonging to the SIVIC. The source is 
detected at 3.6 /xm in only epoch two of the IRAC data, 
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Fig. 6. — SED of source sl 2. T he solid black line shows the best 
fit MARCS model (see Section l3.3l l. The underluminous 3.6 fini flux 
appears to be the cause of the [3.6] — [8] excess apparent in Figure|4] 
and may be due to poor photometry or possibly to absorption from 
acetylene. AKARI finxcs are systematically too bright, which may 
be due to stellar blending. 

and this flux falls well below the apparent stellar contin- 
uum, resulting in what may be an artificially red [3.6] — [8] 
color. 

Despite the under-luminous 3.6 /xm flux and system- 
atically over-luminous AKARI data, the SAGE-SMC 
24 /im measurement of source sl2 does show a very slight 
amount of excess. It is possible that this source may be 
a carbon star belonging to the SMC, and the discrepant 
3.6 ^m flux may indicate variabil ity or possibly molecu - 
lar absorption from acetylene (cf. Ivan Loon et al.l[200l . 
The heliocentric radial velocity for carbon stars on the 
outskirts of the SMC is 149.3 ± 30 km s~^, with a disper- 
sion of Sv ~ 25.2 ± 2.1 km s~^ ([Hatzidimitriou et al.l 
|1997[) . A total of 6% of these carbon stars have velocities 
200 km s^^ < Urad < 250 km s~^. Given that this 
source has a magnitude a nd color very simila r to typical 
carbon stars in the SMC (|Bolatto et al.ll2007l ). it is likely 
not a member of NGC 362, and, as a consequence, no 
signiflcant dusty mass loss occurs below the TRGB in 
NGC 362. 

3.3. The Hertzsprung- Russell Diagram 

In Figure [3 we present a physical Hertzsprung- Russell 
Diagram (HRD) for NGC 362. To determine the stel- 
lar parameters of all stars in the cluster, we flt opti- 
cal, near-IR, and IRAC photometry to the Model At- 
mosp here in Radiative and Convec tive Scheme (marcs) 
code (jGustafsson et al.lll975l l2008f) . The grid of model 
spectra and fl tting technique u s ed he re are described in 
more detafl in lMcDonald et all (|2009f) . To flt the model 
spectra to the data, we give preference to MCPS photom- 
etry, including NOMAD and DENIS photometry only if 
a source is not present in the MCPS catal og. 

A Padova isochrone (|Marigo et al.ll2008l ) is included in 
Figure[7]at 12.5 Gyr and [Fe/H] = - 1.28. This metal- 
licity is sligh tly lower than the [Fe/H] = —1.16 quoted in 
iHarrisI (|1996). which we use to calculate mass-loss rates 
in Section ISTSl We note that a value of [Fe/H] = —1.28 
would result in a 32% increase in the gas-to-dust ratio [iji] 
and therefore a 15% decrease in the dust mass loss. Also 
plotted are isochrones represen ting the SMC and dis - 
tant Galactic stars. As in uj Cen (jMcDonald et al.|[2009l ). 
the NGC 362 upper RGB is cooler than the Padova 



isochrone, which may suggest (dusty or non-dusty) mass 
loss on the RGB. RGB mass loss is not well accounted 
for in Simple Stellar Population models applied to distant 
galaxies, but it may be of critical importance in predict- 
ing effective temperatures and thus colors and bolometric 
corrections (Salaris & Cassisi 2005). 

The brightest/dustiest candidate mass- losing stars in 
the cluster are circled in red in Figure [71 and the mod- 
erately dusty stars are marked with purple diamonds. 
The very red source directly above the galaxy locus in 
Figure [3] (source sll) is underluminous and very cold 
on the HRD, consistent with the assertion that it is a 
background galaxy. Note that two of the moderately 
dusty sources are slightly warmer than the RGB (sources 
s03 and s09), which could indicate that these are p ost- 
AGB stars, ahhough iMcDonald fc van Loonl (120071 ) in- 
dicated photosphere temperatures of « 3900 K for both 
stars, placing them firmly on the RGB. Source s07 is 
the coolest of the mass-losin g stars at T w 3300 K. 
IMcDonald fc van Loonl ()2007[ ) measure a temperature of 
T K, 3950 K, placing s07 among the other strong mass- 
losing stars. All other candidate mass-losing stars have 
typical temperatures and luminosities. 

There is an anomalous source in Figure [7] with T k, 
3700 K and L k, 300 Lq. If this source is at the distance 
of the SMC, its luminosity is L 15 000 Lq , which is 
typical for a carbon star. Source sl2, which we suggested 
may also be an SMC carbon star in Section l3?2l was flt 
to a MARCS model with T w 4800 K and L w 820 Lq 
(see Fig. |6|). However, we note that a lack of photom- 
etry shortward of 0.9 fim for source sl2 resulted in an 
uncertain flt. 

3.4. Distribution of Dusty Stars 

Figure [8] shows the spatial positions of the ten sources 
identifled as candidate dusty mass-losing stars (including 
sOl, a likely cluster nonmember). Within the inner 3.5' 
of the cluster, which includes all of the candidate dusty 
stars, conflrmed cluster members are a mean distance of 
1.2' ± 0.1' to the cluster center. For the dusty stars, 
this is 0.7' ± 0.2'. While the dusty stars are slightly 
more tightly concentrated than the general population, 
the uncertainties easily allow for the possibility that both 
populations have the same central concentration, sup- 
porting the assertion that these stars are indeed cluster 
members. 

Low- mass stars lose most of their mass on the RGB, 
and mass segregation will cause AGB stars to migrate 
to the outer regions of the cluster given that the time 
between the RGB and AGB phases is on t he order of the 
relax ation time scale of NGC 362 (~10^ vr: lFischer et al.l 
119931 ). However, most of the dusty stars in NGC 362 are 
likely RGB stars at the TRGB, and have therefore not 
had enough time for mass loss to alter their orbits, unless 
signiflcant dustless mass loss occurred much earlier. We 
therefore do not expect to find that the dusty stars have 
orbits differing from those of the bulk population. 

3.5. Mass-Loss Rates 

The SEDs of the moderately dusty stars are presented 
in Figure [HI and the SEDs of the mor e extreme stars are 
presented in Figure M AKARI d&ta. Illta et al.ll2007l) are 
plotted as red squares, Spitzer data are plotted as blue 
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Fig. 7. — Physical HRD for NGC 362. As in Figure |4l candidate strong mass-losing stars are circled in red, and candidate moderate 
mass-losing stars (excluding sOl, a likely cluster nonmember) are marked with purple diamonds (see Table |2j- Padova isochrones are also 
plotted in dark blue for NGC 362, light blue (dashed line) for Galactic stars at a representative distance, and magenta (dashed line) for 
the SMC. Note that the RGB is cooler than expected. 



circles, pluses indicate literature optical photometry, and 
best-fit MARCS spectra are plotted with solid lines. None 
of the moderate stars appears strongly variable between 
Spitzer epochs and the AKARI epoch {AKARI data for 
source s09 is unreliable). The extreme stars do show 
potential moderate variability. 

The AKARI data for sources s09 and s07 are very in- 
consistent with the Spitzer data, which may be the re- 
sult of source confusion in AKARI. Without exception, 
AKARI 24 /.tm fluxes are bright compared to Spitzer 
24 /xm photometry. This discrepancy is probably due 
to the poor angular resolution of AKARI compared to 
Spitzer (see Section I2.1|) . resulting in source confusion 
and blending. All but three of t he sources in Table El 
were also detected at 24 ^m by iBarmbv et al.l (|2009( )7 
These 24 pum data are very consistent with the Spitzer 
SAGE-SMC 24 ^m data, typically differing by a tenth 
of a magnitude or less. Source s02 is the only excep- 
tion to this, with a difference of 0.5 m ag between Spitzer 
SAGE-SMC and IBarmbv eTaLl (pOOflh . 

To estimate the mass-lo ss rates, we fit the SE Ds using 
the DUSTY modeling c ode (iNenkova et allll999f ). Follow- 
ing iMcDonald et alj (|2009f) . the SEDs are fit by mod- 
els usi ng dust comp osed of either amorphous carbon 
(AMC. iHanned [l988n or silicates, with the latter com - 
posed of 65% astronomical silicates (jDraine 
15% compact AI2O3 (optical constants from the Jena 
database ^), and 10% ea ch of glassy and crystalline sili- 
cates (Ja ger et al.lll99^ . These two dust compositions 
result in good fits, but the dust may consist of other sili- 
cate or carbon species or different proportions of species. 
Carbonaceous dust yields relatively blue [8] — [24] colors, 

* http:/ /www. astro. uni-jena.de/Users/database/entry.html 



whereas oxygen-rich dust always produces a relatively 
large 24 ^m excess compared to 8 /xm, although without 
AKARI ddXd. to fill the gap between 8 and 24 /xm left by 
Spitzer^ it is difficult to definitively distinguish between 
different du st compositions. 

Following IMcDonald et al] (|2009f ). the resulting mass- 
loss rates from DUSTY were converted to real mass-loss 
rates using the following prescription, which implicitly 
assumes a scaling of the wind speed based upon a dust- 
driven wind formalism: 



-^DUSTY 

200 



104 i 



3/4 



200/ 



-1/2 



(2) 



where ip is the gas-to-dust ratio (-0 = 2891, assum- 
ing [Fe/H] = -1.16 and = 200), L is the stellar 
luminosity determined from the MARCS best fit, and 
P(j is the bulk grain density. For silicates, we assume 
Pd = 3 g cm^^, and for AMC grains, we assume pd = 
2.5 g cm"'^. Note that a decrease in metallicity to [Fe/H] 

— —1.28 would result in a 15% decrease in Mdust- A sum- 
mary of input model parameters and output mass-loss 
rates is presented in Tabled 

The moderately dusty stars are all fit well by sili- 
cate dust. In addition, source s09 is fit equally well by 
AMC dust due to a lack of rehable AKARI data. The 
mass-loss rate for s09 from the best-fit AMC model is 
a factor of 1.7 less than the rate from the best-fit sili- 
cate model. Together, the four moderately dusty clus- 
ter members lose Mdust = 8.3(±0.6) x 10"1°Mq yr'^ 
where the errors come solely from the dust composi- 
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Fig. 8. — Distribution of candidate mass-losing stars. Most dusty 
stars are located within the inner 1' of the cluster center. Red 
circles mark sources with strong mass loss, and yellow diamonds 
mark sources experiencing moderate mass loss. The open blue 
triangle (source sOl) is a likely cluster non-member. Source sll, a 
background galaxy, falls well outside the field of view shown here. 

tion used for source s09. This corresponds to Afgas = 
2.4(±0.2) X IO-^Mq yr-i. 

It is not surprising to find silicate dust in the o u tflows 
of t hese four mod erate stars. iHofner fc Andersenl ()2007D 
and iHofneH ()2008il suggest that winds could be driven by 
large (~ 1 iim) Fe-free silicate grains, or that silicates can 
easily form in winds that are driven by carbon grains. 

3.5.1. The most extreme dust producers 

Five sources are classified as strong mass-losing stars 
(s02, s05, s06, s07, and s08). We do not attempt to fit 
a DUSTY model to source sOS due to inconsistent Spitzer 
and optical photometry and a lack of AKARI photome- 
try (Fig. Illjl. The remaining four stars do appear mod- 
erately variable in the IR between Spitzer and AKARI 
epochs (Fig. [TU)) . with s05 showing the largest ampli- 
tude followed by s02 and s06, which are the known vari- 
able stars V2 and V16, respectively. The total dust- 
traced mass-loss rate from the four extreme stars is 
Mgas = 6.2(±0.5) X IQ-^Mq yr-\ with the uncertainty 
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s02 


500 


AMC 


5.9 X 10- 


10 


1.7 X 10" 


-6 




s03 


1200 


silicate 


2.5 X 10- 


10 


7.1 X 10- 


-7 




s04 


1200 


silicate 


2.1 X 10- 


10 


6.1 X 10" 


-7 




sOS'' 


1000 


silicate 


3.2 X 10- 


10 


9.3 X 10- 


-7 




s06 


600 


AMC 


7.0 X 10- 


10 


2.0 X 10" 


-6 




s07 


800 


silicate 


4.6 X 10- 


10 


1.3 X 10" 


-6 




sog'^ 


1200 


silicate 


2.6 X 10- 


10 


7.4 X 10" 


-7 




slO 


1200 


silicate 


1.7 X 10- 


10 


4.8 X 10" 


-7 



TABLE 3 

Mass-loss Rates 



Note. — Sources s02 and s06 are fit equally well to a 
DUSTY model composed of a combination of silicates and 
AMC dust, resulting in a nearly identical mass-los s rate 
to that derived from the AMC models. See Section [3.5.2l 
for a description of the uncertainties in the mass-loss 
rates quoted here. 

^ The silicate model fits well to the Spitzer daXa for s05. 
AKARI da.ta, for s05 is better fit by an AMC model with 
Tinner = 800 K, Mjust = 5.6 X lO'iO Mq yr-i, and 
Afgas = 1.6 X 10-6 Mq yr-i. 

^ Source s09 is fit equally well to a model with AMC 
dust and Tinner = 800 K. The resulting mass-loss rates 
are a factor of 1.7 less than the rates derived from the 
silicate model. 

depending on the dust compositions chosen. This is more 
than twice the contribution of the four more moderate 
stars. 

The dust composition is difficult to estimate for the 
four extreme sources. Source s05 is particularly odd. 
The AKARI and Spitzer data appear to show two dif- 
ferent dust compositions. AMC dust fits best to the 
AKARI data, but is simultaneously ruled out by the 
Spitzer data due to a lack of e xcess between 3.6 and 
8 ^m. iMcDonald k, van LoonI (|2007l ) found that this 
source is actually two stars (o05a/o05b) with similar 
temperatures and luminosities (T « 4100 K, Lqosq = 
2103 ± 189 LQ,Lo05b = 2392 ± 213 Lq), which are to- 
tally unresolved in both Spitzer and AKARI data. Both 
stars appear oxygen-rich and show blue-shifted Ha ab- 
sorption cores that indicate an outflow of material. It 
is unclear how their small separation (0.2") has affected 
the IR photometry. It is possible that the low resolu- 
tion of the AKARI images has resulted in blending of 
sources in addition to o05a and o05b, causing inflated 
AKARI fluxes. Another possibility is that one or both 
of these sources is variable, producing different types of 
dust in varying proportions during different phases of the 
pulsation cycle. A third possibility is that neither star 
has a real IR excess, and the apparent excess is caused 
by inaccurate flux extraction due to stellar blending in 
this crowded region of the cluster. However, this third 
possibility seems unlikely since s05 also shows a strong 
24 iim excess (Fig.[TU]), suggesting that at least one of the 
stars is producing dust. Nevertheless, the Ha line clearly 
shows that both stars (o05a/o05b) are losing mass, al- 
though it may not be as clear that the mass loss in either 
star is accompanied by dust production. 

For source s07, AKARI photometry is available only 
at 24 /.tm, and this flux is more than twice the flux mea- 
sured with MIPS at 24 ^m. If the AKARI 24 /zm flux 
is accurate, it can only be explained by an extremely 
large quantity of circumstellar dust. This, along with the 
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Fig. 9. — SEDs of moderately dusty stars at the TRGB in NGC 362. Wavelength error bars correspond to the filter width. AKARI data 
are plotted as red squares and Spitzer data are plotted as blue circles. MARCS models are solid black lines. IRAC epoch 1 is plotted with 
open blue circles, and IRAC epoch 2 is plotted with closed blue circles. All four stars are well fit by a DUSTY silicate model (magenta, 
dotted line), but s09 is fit equally well by a DUSTY amorphous carbon (AMC) model (green, dashed line). 



source's location near several other bright 24 /xm sources, 
suggests that the AKARI 24 /im flux is overestimated. 
We therefore exclude the AKARI 24 /.tm point from the 
DUSTY fit and find that silicate dust provides the best fit 
due to the lack of strong excess in the IRAC data. 

The two remaining extreme sources, s02 (V2) and s06 
(VI6/0OI), are the most extreme stars in the cluster and 
exhibit the strongest mass-loss rates. Neither star is fit 
well by silicates alone, but both are well fit by AMC 
models. This is surprising, give n that the optical spec- 
truni o f s06 indicates C/0 < 1 ([McDonald fc van LoonI 
I2007| l. I Smith et all (|1999l ) did not measure the carbon 
abundance in star s02, so its nature as carbon or oxygen- 
rich remains unknown. C arbon stars are knowii to reside 
in GCs, notably in lo Cen (|van Loon et al.ir2007[ ). Models 
that include a combination of 50% AMC dust and 50% 
silicate dust also fit well to both s06 and s02, but we 
note that to fit this model to s02, we must exclude the 



^Tfyli?/ photometry, which shows no silicate features. 

A mid-IR spectrum of s02 was obtained by Y. Ita et al. 
(in preparation) using AKARIlYiC slit-less spectroscopy. 
Pre-analysis of the spectrum suggests a complete lack of 
silicates and the possible presence CO and water absorp- 
tion near 3 and 5 /^m, which may explain the underlumi- 
nous photometry near these wavelengths. While 

we cannot definitively declare a dust type in either s02 
or s06, AMC dust must play some role if we are to rec- 
oncile the strong photometric excess in the IRAC bands 
and the excess at 24 /im. 

3.5.2. Uncertatnttes in the Mass-Loss Rates 

The uncertainties quoted for the total mass-loss rates 
of the candidate mass-losing stars reflect a range of val- 
ues that depends on the dust composition chosen. The 
uncertainties in the mass-loss rates for individual stars 
is not reflected in this value. In uj Cen, uncertainties of 
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Fig. 10. — Same as Figure[9l but for the most extreme stars in NGC 362. The amorphous carbon (AMC) + siUcate DUSTY models shown 
for s02 and s06 are a combination of 50% AMC and 50% silicates. 
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Fig. 11. — SED for source s08. We did not attempt to fit a 
DUSTY model to this SED due to inconsistent Spitzer and optical 
photometry and a lack of AKARI photometry. This source may 
be affected by stellar blending or it may be extremely variable. 

52 - 80% were derived for the mass-loss rates of individ- 



ual stars, based on internal errors from DUSTY and on 
several assumptions. These assumptions include (1) pul- 
sation causes only negligible variations in the mid-IR, (2) 
the gas-to-dust ratio scales with metallicity, (3) the wind 
velocity scales as Wwind cx; 

7 ,1/4 1/2^ and ( 4) du st-to- 
gas coupling is efficient. See lMcDonald et all ()2009D for a 
more detailed discussion of mass-loss rate uncertainties. 

We note that the prescription we use to estimate the 
mass-loss rates assumes that the wind speed scales as 
"^wind L^/'*(-!/)pd)~^^^, resulting in velocities ranging 
from Wwind — 0.5 — 1.3 km s^^ for mass-loss candidates 
in NGC 362. These small velocities provide the largest 
uncertainty to our mass-loss rates, as they are on the 
order of the turbu lence in the wind ('-~^1 — 2 km s^^; 
iSchdier et al.l[200^ . It is not clear whether very slow 
winds can be sustained, but it is possible that wind 
speeds can be increased to <10 km s~^ due to the 
effect s of rotation and puls a tion in the star. More- 
over, iMcDonald fc van LoonI ()2007l ) estimated that the 
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wind velocities in several NGC 362 stars range from 
«6 — 17 km s~^. If these velocities are good approxi- 
mations of the wind speed in the dust-producing zone, 
then the mass-loss rates of each star will increase by up 
to a factor 6. 

4. DISCUSSION 

4.1. Dust and Mass Loss in V2 and V16 

Star s06 (V16) appears to show carbon dust despite its 
nature as an oxygen-rich star. Star s02 (V2) is similar, al- 
though its nature as carbon- or oxygen- rich is unknown. 
Thes e stars are akin to V42 in oj Cen (cf. IMcDonald et al.l 
12001 . an oxygen-rich star which once may have shown a 
strong 10-/im silicate dust feature, but more recently ap- 
pears to show predominantly amorphous carbon dust de- 
spite its oxygen-rich atmosphere. While uj Cen V42 and 
NGC 362 V16 both appear to be currently dominated 
by carbonaceous dust, they may still contain some (pos- 
sibly variable) amounts of silicates. 47Tuc is home to 
two dusty star s (VI and V18) th a t show variable silicat e 
emission (van Looii et al.ll2006al : iLebzelter et al.ll2006D . 
iHofner k, AndersenI (|2007l ) suggest a possible scenario 
where non-local thermodynamic equilibrium conditions 
allow the periodic production of different grain composi- 
tions at different epochs in the pulsation cycle. It may 
also be the case that grain size varies enough that large 
grains periodically suppress the silicate features. 

4.2. Global Mass Loss 

The cluster stars that exhibit IR excesses are cur- 
rently returning Mgas = S.eig^ x 1O"'^M0 yr"^ of 
ga s to the ICM, us i ng th e same fractional errors used 
in IMcDonald et all (|2009f) . To find the total cluster 
mass loss (dusty and dustless), we estimate the mass- 
loss rates for no n-dusty stars usi ng the method de- 
scribed in Schrod er fc Cuntd (120050. which is based on 
the stellar temperature and luminosity. We find that 
chromospherically-driven (dustless) mass loss is respon- 
sible for 2 X lO^^Af0 yr^^, giving a total cluster mass 
loss of A'/totai ~ l-ltal X lO^^M© yr-^ Dust produc- 
tion therefore accompanies 80% of the current mass loss 
in NGC 362. 

The global mass-loss rate for NGC 362 is virtually 
identical to that of u ;Cen (w 1.2 x 10"^ Mq yr^^; 
IMcDonald et al.ll2009f) . despite NGC 362 being an or- 
der of magnitude less massive than the latter. 47 Tuc, 
which is similar in mass to oj Cen, appears to have a 
similar global mass-loss rate, given that its four long- 
period variables h ave mass-loss rates ranging from 10~^ 
to 10-^Mq yr-i (Origl ia et aIll2007D . In aU three clus- 
ters, only a handful of stars are dominating the global 
mass-loss rate. This supports the likelihood that global 
dust production and dust-associated mass loss in clusters 
are stochastic in nature, driven by episodic dust pro- 
duction in the individual stars (possibly due to varying 
stellar conditions at different pulsation phases) and the 
rarity of the stars themselves, lo Cen may currently be 
experiencing a relatively quiescent period of dust pro- 
duction, given that only two stars (V6 and V42) are cur- 
rently forming significant amounts of dust. M15, on the 
other hand, appears to have recently experienced a pe- 
riod of strong dust production, assuming that the ICM 
dust present there was formed in stellar winds rather 



than by a stellar collision (|Rasio fc ShapirolfTMl) or one 
or more diffusing planetary nebulae. NGC 362 has a 
moderate total mass (3.78 x 10^ -^0) and up to four 
stars producing significant amounts of dust. NGC 362 
may therefore be representative of a more intermediate 
phase of dust production in the typical GC. 

Any effects that metallicity has on dust production 
appear to be overshadowed by the apparent episodic 
nature of dust production. Stars with strong IR ex- 
cess are discovered in clusters with metallicities ranging 
from [Fe/H] = -2.4 for M 15 ("e.g.. [Bover et al. 2006) to 
[Fe/H] = -0.7 for 47 Tuc (|van Loon et al.ll2006aD . with 
no clear trend favoring higher metallicities. 

Dust production is dominated by the most luminous 
stars in the cluster. However, several dustZess stars 
have temperatures and luminosities similar to the dust- 
producing stars. The Ha profiles of at least a handful 
of these dustless stars suggest that they are losing mass 
without any associated dust production. The only ob- 
served difference known to exist between the dusty and 
non-dusty stars is a lack of a red Ha line emission wing in 
the non-dusty stars, perhaps suggesting that material is 
rest ricted to smaller stellar radii than in the dusty stars 
(cf. IMcDonald fc van Lom]|[200l . 

4.3. The Intracluster Medium 

If we make the (possibly naive) assumption that the 
rate of gas return has remained relatively constant since 
NGC 362 last plunged through the Galactic plane (3x lO"^ 
yr ago), then we might expect 210 - 540 Mq of gas and 
0.05 - 0.15 Mq of dust to have gathered in the ICM, 
especially considering that the cluster escape velocity is 
much larger than typical wind speeds derived with dusty 
(wcsc,o — 46.7 km s^^, Wwind — 0.5 — 1.3 km s"-'^, assum- 
ing Uwind oc L^/^(-(/;pd)~^^^)- Rotation and magnetic ac- 
tivity and pulsation shocks can incr ease wind ve locity 
by up to an order of magnitude (cf. lBowen|[l988D . but 
this increase is still not sufficient for material to escape 
the cluster. Within the half-ligh t radius, upper lim its of 
only < 6.3 x IO'^Mq of dust (Barmby et al. '2009) and 

< 1.8 Mq of ionized hydrogen (Hesscr fc Shawl 19715) 
have been determined in the ICM of NGC 362. 

Assun iing that 5.5% of t he ICM gas is ion ized, as in 
47 Tuc (jSmith et all [l990l : iFreire et al.l [200lf) . we find 
that there is < 40 Mq of neutral gas in the ICM of 
NGC 362. Based on the total (dust + dustless) clus- 
ter mass-loss rate (1.1 x 10~^ Mq yr~^), we can expect 
neutral gas to be cleared from the cluster on timescales 

< 3.6 X 10® yr, more than an order of magnitude shorter 
than the time since the last Galactic plane-crossing. 
The dust provides an even shorter timescale of ICM 
removal or destruction, namely < 2 x 10^ yr, given a 
dust mass- loss rate of 3 x 10~^ Mq yT~^ and a dust 
upper limit of < 6.3 x lO~^Af0. To escape the half- 
mass radius (0.81') within 20 000 yr, the dust would 
have to travel > 80 km s~^. This velocity is much 
higher tha n the likely w ind velocities (t'wind 5. 10 km s~^, 
IMcDona l d fc van Loonj [2007; Me szaros et al.ll2009[ ). re- 
quiring that (a) either the dust is dissociated on short 
timescales within the cluster or that (b) dust is rapidly 
accelerated out of the cluster. 

Based on the space velocity of NGC 362, (140 km s~^, 
cf. lOdenkirchen et al.. ,1997i) . ram-pressure from hot 
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Galactic Halo gas could clear the cluster ICM in «2.8 x 
10^ yr. A stellar collision has enough kinetic energy to 
clear the ICM, and such collision s occur on timescale s 
of «2.5 X 10*5 yr in NGC 362 (cf. iBarmbv et aTI[2009f ). 
If ICM dust were destroyed within the cluster on the re- 
quired shorter timescales, then ram-pressure stripping or 
a stellar collision may be responsible for the lack of gas 
and the byproducts of dust destruction in NGC 362. 

Although mass loss estimates predict the presence 
of ~10^ Mq of ICM material in most GCs, nothing 
approaching this amount of mass has ever been ob- 
serve d. a; Cen has an upper hmit of < 10~* Mq of 
dust (iBover et al. l l2008D and < 2 8 of neutral gas 
(|Smith et al.lll990D . IBarmbv et all ((200a) finds ICM dust 
mass upper limits of < 10"'* AIq in eight GCs, and 
iMatsunaga ct al. ( 2008) finds upper limits of < 10~^ Mq 
of dust in 12 GCs. van L oon et al.l (|2009f) determined 
3(7 upper limits between 6 and 51 Mq of neutral hydro- 
gen within the tidal volumes of four GCs. Much of the 
ICM may be ionized, as demonstrated in 4 7Tuc. where 
up to 0.1 Mq of ionized gas was discovered (jFreire et al.l 
1200 It) . M15 shows the strongest evidence for an ICM, 
with firm detections of 9 x 10~^ Mq of d ust (IBover et al 



20061) and 0.3 Mq of neutral hydrogen (jvan Loon et al 



2006bD . The combination of episodic dust production 



and a stochastic ICM removal process (such as a stel- 
lar collision) could explain the lack of dusty ICM ma- 
terial in GCs and its unusual presence in M15. If 
ram-pressure stripping from hot Halo gas is the domi- 
nant removal mechanism, the intense crowding in M15 
(log(rt dai/llcore) = 2.5, Compared to 1.9 for NGC 362; 
IHarrisl[l99l could cause the stripping of material in the 
center of the cluster to be less efficient, with bow-shocks 
around outer stars effectively shielding the cluster's ICM. 

It has become clear that low-mass, low-metallicity stars 
successfully form dust and lose mass. Therefore, some 
mechanism must remove gas and dust from GCs. This 
material may ultimately reside in the Galactic Halo, 
fall back onto cluster st ars, or find its way back into 
the Milky Way disk (e g..lEyans et al.l l2003l: feover et all 
l2?)06t I van Loon et al. I I2006H I2009f r" In any case, low- 



mass, low-metallicity stars like those in globular clus- 
ters must collectively and continuously contribute a large 
amount of recycled and enriched material to the ISM, 
helping to drive galaxy evolution. 

5. SUMMARY OF RESULTS AND CONCLUSIONS 

We present an analysis of dusty mass loss in the 
Galactic globular cluster NGC 362 carried out with 
data serendipitously obtained during SpitzerS AGE-SMC 
observations of the Small Magellanic Cloud. Spectral 
energy distribution modeling of all cluster stars pro- 
vided stellar parameters, allowing the construction of 
a physical Hertzsprung-Russell Diagram. The Red Gi- 
ant Branch (RGB) is slightly cooler than the Padova 
isochrone, indicating significant mass loss on the RGB. 

Significant infrared (IR) excess exists only at and 
above the tip of the RGB. The four brightest stars at 
24 /im exhibit the strongest IR excess. Four additional 
stars near the tip of the RGB show moderate IR ex- 
cess. Modeling with the dusty code indicates that these 
eight stars contribute 8.6I3 4 x IO'^Mq yr~* of gas 

and S.Otl^ x 10"^Mq yr"* of dust to the intracluster 
medium. The two most extreme stars, variables V2 and 
V16, provide up to 45% of the total (dust accompanied) 
mass loss. 

The eight dustiest stars show evidence of silicate dust, 
but, surprisingly, three of the four more extreme stars re- 
quire some amount of amorphous carbon dust to explain 
their mid-IR excesses. Strong mass loss therefore appears 
to correlate with a larger contribution from carbon dust, 
which may suggest that non-equilibrium conditions are 
common in such stars. 
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